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Abstract 

We provide a systematic construction of string junctions in curved backgrounds which are 
relevant in computing, within the gauge/gravity correspondence, the interaction energy 
of heavy dyons, notably of quark-monopole pairs, in strongly coupled SYM theories. 
To isolate the configurations of physical interest we examine their stability under small 
fluctuations and prove several general statements. We present all details, in several 
examples, involving non-extremal and multicenter D3-brane backgrounds as well as the 
Rindler space. We show that a string junction could be perturbatively stable even in 
branches that are not energetically the most favorable ones. We present a mechanical 
analog of this phenomenon. 
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1 Introduction 



String junctions [1, 2] arc of particular interest in string theory for several reasons: Due to 
their BPS nature they can be used to build supersymmetric string networks [3, 4, 5, 6, 7] 
and they are useful for the description of BPS states in SYM theories [8] and for appli- 
cations to gauge symmetry enhancement in string theory [9]. The scattering of string 
modes at string junctions in flat spacetime has been analyzed in [10]. Particularly, impor- 
tant for the purposes in the present paper is the fact that since string junctions connect 
strings of different type they can be used, within the gauge/gravity correspondence [11], 
to compute the interaction energy of dyons and in particular of heavy quark-monopole 
pairs [12]. 

The purpose of the present paper is to formulate and construct the string junctions ap- 
propriate for computing the interaction potentials of heavy dyons at strong coupling, 
within the gauge/gravity correspondence. This is done for a large class of curved back- 
grounds away from the conformal point and with reduced or no supersymmetry at all. 
Having constructed string junction solutions does not imply that the corresponding dyon 
interaction potentials are physical. One has to at least investigate perturbative stability 
which will appropriately restrict the parametric space in which the solutions are physi- 
cally relevant. Such investigations were exhaustively performed [13, 14] for single string 
solutions useful in computing the heavy quark- ant iquark potential. Given these works 
and in comparison to them we will find that the stability analysis of junctions leads 
to expected as well as unexpected results, the latter being counterintuitive to stability 
arguments based on energy considerations. 

The organization of this paper is as follows: In section 2 we formulate string junctions 
in a quite general class of curved backgrounds. We derive general formulae from which 
the interaction energy of a dyon pair is computed. In section 3 we consider in great 
detail the stability analysis under small fluctuations of these string junctions aiming at 
discovering the physically relevant regions. We pay particular attention to the precise 
formulation of the boundary and the matching conditions at the junction point and derive 
general statements that isolate the boundaries of stability for all types of perturbative 
fluctuations. In section 4, we present several examples of string junctions using D3-black 
branes and multicenter D3-brane solutions which within the AdS/CFT correspondence 
are relevant for jV" = 4 SYM at flnite temperature and at generic points in the Coulomb 
branch of the theory, respectively. In addition, we present the details of string junctions 
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in Rindlcr space, given its relevance in a variety of black hole backgrounds near the 
horizon. In section 5, we apply the outcome of the general stability analysis of section 3 
to the examples of section 4. We present our conclusions in section 5. In the appendix 
work out the stability of a classical mechanical system which also exhibits an unexpected 
perturbative stability similar to the one we found with string junctions. 

2 The classical solutions 

In this section we develop and present the general setup of the AdS / CFT calculation of 
the static potential of a heavy static dyon-dyon pair with general NS and RR charges. 
These computations involve three strings and therefore we will be based to the well-known 
results for the potential of heavy quark- antiquark pairs of [15] for the conformal case, 
as extended in [16] for general backgrounds, and to the computation of the interaction 
energy of a heavy quark-monopole pair in the conformal case [12]. 

We consider a general diagonal metric of Lorentzian signature of the form 

Here, y denotes the (cychc) coordinate along which the spatial side of the Wilson loop 
extends, u denotes the radial direction playing the role of an energy scale in the dual 
gauge theory and extending from the UV at it — > oo down to the IR at some minimum 
value Mmin determined by the geometry, x stands for a generic cyclic coordinate, 9 stands 
for a generic coordinate on which the metric components may depend on and the omitted 
terms involve coordinates that fall into one of the two latter classes without mixing terms. 
It is convenient to introduce the functions 

g{u,9) = —GttGuu , fy{u,9) = —GuGyy , fx{u,9) = —GuGxx , 
fe{u,9) = -GttGee, h{u,9) ^ GyyG^u ■ (2.2) 

If the conformal limit can be taken (this is the case in all examples in the present paper 
except the one involving the Rindler space) we can approximate the metric by that for 
AdS^ X with radii (in string units) R — {AngsNy/"^, with Qs being the string coupling. 
In this limit we have the leading order expressions 

gr~/i~l, fy'^u^/R'^ , fe^u^, as u ^ oo . (2.3) 
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A string junction consists of three co-planar strings joined at a point as depicted in Fig. 
1. The string labelled as 1 has {p,q) NS and RR charges, the second labelled as 2 has 
{p',q') charges, whereas the third straight string labelled by 3 has {m,n) = {p + p',q + q') 
due to the charge conservation at the junction point. From a microscopic point of view 
the strings that interact are the fundamental F- and D-string with charges (1,0) and 
(0, 1), respectively. Any other charge combination should be achievable by performing in 
this pair an 5'L(2, R) transformation which gives the condition p q' — q p' — ±1 [9]. The 
first two strings end at the boundary of the space-time at u — oo, whereas the straight 
third string ends at the minimum value itmin of the space-time allowed by the geometry. 
All three strings meet at the junction point at u = uq. Within the framework of the 
AdS/CFT correspondence, the interaction potential energy of the (p, q) dyon is given by 

g-iBT ^ (^w{C)) = e'^"^"^^ , (2.4) 

where 

Sp,q[C] = J rfrda^^det^ , g^p = G^^dax'^d^x" , (2.5) 

is the Nambu-Goto action for a string propagating in the dual supergravity background 
whose endpoints trace a rectangular contour along one temporal and one space direction, 
Tp^q — _|_ q'^/g'^ and similar expressions hold for the strings 2 and 3. We will not 
consider a contribution from the Wess-Zumino terms either because these terms do not 
exist all or have vanishing contribution within our ansatz for solutions below. 

We next fix reparametrization invariance for each string by choosing 

t — T , u — a , (2.6) 

we assume translational invariance along t, and we consider the embedding 

y = y{u) , x = , 9 — 9o — const. , rest = const. , (2.7) 

supplemented by the boundary condition 

u{-Li)=u{L2)=oo , (2.8) 

appropriate for a {p,q) dyon placed at y — —Li and a ip',q') dyon placed at y — L2. 
The string with charges (m, n) extending from the junction point to u — Umin is straight 
as y = const, is always a solution of the equations of motion (see (2.11) below). In the 
ansatz (2.7), the constant value ^0 of the non-cyclic coordinate 9 must be consistent with 
the corresponding equation of motion. This requires that [13] 

deg{u,9)\e=e, = d0fy{u,9)\e=eo = . (2.9) 
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Figure 1: Junction configuration: We denote by Fj, i = 1,2,3 ( Fpg etc in the text) the 
forces exerted at each string at the junction point uq. The turning points Ui, i = 1,2 for 
the strings 1 and 2, are not depicted and he at the right of Mq. We always have Mmin < Uq, 
but the Mj's may or may not be larger than Wmin- 

For the ansatz given above, the Nambu-Goto action reads 

T T 

n n 



^■^ ~ 27r 



J du^g{u) + fy{u)y'^ , (2.10) 



where T denotes the temporal extent of the Wilson loop, the prime denotes a deriva- 
tive with respect to u while g{u) = g{u,9o) and fy{u) = fy{u,9o) are the functions in 
(2.2) evaluated at the chosen constant value 6q of 6. Similar actions are also considered 
for the (p', q') string, as well as for the straight (m, n) string. Independence of the La- 
grangian density from y implies that the associated momentum is conserved, leading to 
the equation 

^ y'^,^^^M^ (2.11) 



where Ui are the values of u at the turning point for each string, fyi = fyiui) , fyo = fy{uo) 
and yci is the classical solution with the two signs corresponding to the lower (string 1) 
and upper (string 2). The symbol Fj (not to be confused with the forces in Fig. 1) stands 
for the function 

Fi^j^, i = l,2. (2.12) 

Jy Jyi 

Note that, for the junction point uq we have that uq ^ Mmin. For the turning points we 
have that Ui ^ uq, but they are not necessarily larger than Mmin, since the strings 1 and 
2 are not actually extending to their turning points. Integrating (2.11), we express the 
separation length as 



oo fT^ poo pn^ 

r\ , „l/2 / , V-f'2 

Wo jy •'"0 Jy 



L, + L, = fli^ / du^ + 4f / du^ . (2.13) 



Finally, inserting the solution for y^j into (2.10) and subtracting the divergent self-energy 
contribution of disconnected worldsheets, we write the interaction energy of the junction 
as 

noo poo 

£i= du^/Fi- du^ , 1=1,2. (2.14) 



Uo 



Ideally, one would hke to evaluate the integrals (2.13) and (2.14) exactly, solve (2.13) 
for Uq and insert into (2.14) to obtain an expression for the energy E in terms of the 
separation length L and Ui, U2. Then, minimizing this expression with respect to the Wj's 
we will obtain the expression of the E in terms of the separation length L, which is then 
identified with the interaction energy of the heavy dyon pair. 

However, in practice this cannot be done exactly, except for the conformal case, and Eqs. 
(2.13) and (2.14) are to be regarded as parametric equations for L and E with parameters 
Mo, Ml, U2. A much easier way to find mi, U2 is to impose that the net force at the string 
junction is zero [2, 4, 10], an approach followed in the present context for the conformal 
case in [12]. The infinitesimal lengths squared along each string are 

de^ = {Gyyy'2 + Guu)du^ = -^Fidu^ , i = 1, 2 . (2.15) 

Hence from the action (2.10) we find that the tensions of the strings at the junction point 

IP 

u = Uo are \/ —Gtt and similarly for the other two string that meet at the junction. 
27r 

The angles between the string and the u axis at the junction point are computed from 

\fG^yd.y ^ 



At M = Mo : tan di 



VG^du 




sinOi = , i = 1,2 , (2.16) 



where we used (2.11) and took into account that displacements along perpendicular axes 
arc measured by a curved metric. The above explicitly shows that the angles depend on 
the parameters Uq and Mj. The forces exerted to each of the strings on the u—y plane are 

Fp,q = -^V-Gu (-cos^i,sin^i) , 

Fp',,' = -^^^ (cos^2,sin^2) , (2.17) 

Fm,n = -^^V-Gtt (1, 0) , 
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the first entry being the u and the second the y component. Demanding that the total 
force be zero one finds that the angles at the equilibrium point are given by 



cos Q\ — 



-m,n ' p,q p ,q 
m,n-'-p,q 



rp2 _|_ rp2 rp2 

COS 02 = . (2.18) 

■^-^ m,n-'-p',q' 

Prom these expressions the angles are determined in terms of the NS and RR charges of 
the strings and then from (2.16) one may express Ui in terms of Uq and the string charges. 
Recalling that Uq is determined by L, one ends with the interaction energy being solely 
a function of L and of the strings' charges as fixed parameters. 

We can verify that the approach of fixing the parameters Ui by utilizing the zero force 
condition at the junction point gives rise to a minimal energy. Since we can not solve 
(2.13) for Uq in general, we shall recall that we have an implicit expression of uq in terms 
of L,ui,U2, by inverting (2.13). Since we choose to minimize the energy with respect to 
Ui, the physical length L does not depend on them and we have 

dL 

^ = 0, i=l,2, (2.19) 

which can be used together with (2.16) to express the derivative of uq with respect to Ui 
as 



^""^ - ^ I du-y^^J^ , i = 1, 2 . (2.20) 



dui tan 6'i + tan 6^2 2 ]] gofyi Juo ify-fyiY''^' 

We have not used the no force condition to derive the above in the sense that the angles 
9i are not constant except for the equilibrium point. Next, using (2.20) and the zero 
force condition one can check, after careful algebraic manipulations, that 

dE 

— = 0, ^ = l,2, (2.21) 

OUi 

as advertized. Our proof is completely general within our ansatz (2.1) and no- where we 
did need to explicitly evaluate the various integrals by resorting to specific examples.^ 
Had we done so, in the cases where this is possible, we would have ended up with 
complicated expressions and our proof would have been equivalent to proving certain 
identities involving, in general, special functions (for instance, for the conformal case see 
the app. of [12]). 



^The result is also valid even when in the upper limit of integration infinite is replaced by a finite 
value, as in the example of the Rindler space below. 
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Let us mention that the expression for the interaction energy E in terms of length L for 
the configuration of {p, q) and (g, p) dyons is invariant under the S-duahty transformation 
Qs ^1 Qs-i as it should be. However, this is not true for L and E separately as functions 
of the auxiliary parameter uq. 



2.1 On the quark- monopole interaction 

The most important case arises for the interaction potential of a quark with a monopole 
with charges (1, 0) and (0, 1), respectively. In that case it is easily seen from (2.18) that 
9i -\- 62 — 'K I'i- It will be useful for later purposes to consider the limit of — > 0. In 
that case the tension of the (0, 1) string corresponding to the monopole becomes very 
large, so that we expected that it stays almost straight, whereas in order for the forces 
to balance, the (1,0) string should hit the straight string almost perpendicularly. This 
is indeed reflected in the expansions 

01 = \-9s + \9l + 0{gl) , 

02 = Qs-lg's+Oig!) . (2.22) 

Of course an equivalent expansion exists for gs ^ 00, with 9i and 62 interchanged and 
gs — > l/s's- In the small string couphng limit we easily infer from (2.16) and (2.22) that 
the turning points of the strings 1 and 2 are 

ui^uo + 0{gl) , U2 = liroot + 0{gt) , (2.23) 

where Wroot is the largest root of the equation fy2 = 0. In the examples we present in detail 
below, itroot — i^min, except for one case in which it has a smaller value. Therefore the 
expressions for the separation length (2.13) and the energy (2.14) for the quark- monopole 
interaction, become 

,1/2 / , 



and 



fyo I + 0{gs) (2.24) 

jy 



uo Jy 



E = ^[ duy/Fo - / du^ + 0{gs) , (2.25) 

where we have defined Fq as in (2.12) with fyi replaced by fyo = fy{uo). Thanks to 
the ^-duality invariance of the general results (2.13) and (2.14) identical expressions to 
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(2.24) and (2.25) hold in the large string coupling expansion, where in the correction 
terms we just replace Qs 1/fi's- Hence, in the limit of small and large string coupling, 
the energy and the length assume half of their corresponding values for the single string 
for the quark- antiquary system. We emphasize that this fact does note imply that the 
spectrum of small fluctuations around the each classical configuration will be the same. 
Actually we will explicitly demonstrate that the opposite is true, which implies that no 
matter how stiff two of the strings (the straight and one of the strings 1 or 2) become, 
its fluctuations do not decouple and affect those of the third string. 

3 Stability analysis 

We now turn to the stability analysis of these conflgurations, aiming at isolating in para- 
metric space the physically stable regions and pointing out the similarities and differences 
with the similar analysis for string configurations dual to the quark- antiquark potential. 
The small fiuctuations about the classical solutions we will discuss, fall into three types: 
(i) "transverse" fluctuations, referring to cyclic coordinates transverse to the dyon-dyon 
axis such as x, (ii) "longitudinal" fluctuations, referring to the cyclic coordinate y along 
the dyon-dyon axis, and (iii) "angular" fluctuations, referring to the special non-cyclic 
coordinate 6. In [13] general results for single string fluctuations of all of the above 
types for the class of backgrounds with metric (2.1) were derived and used to analyze 
the stability of heavy quark- antiquark potentials. These are certainly relevant for the 
investigations of the present paper, so that we review them below, but one has to be 
careful by paying particular attention to the matching conditions at the junction point of 
the three strings which affect the parametric space where the stability occurs, in crucial 
ways. 

3.1 Small fluctuations 

We parametrize the fluctuations about the equilibrium conflguration, by perturbing the 
embedding according to 

Xi = 5xi{t, u) , yi = yci,iiu) + 6yi{t, u) , 9i = 9o + d9i{t, u) , (3.1) 

where the index i refers to the three strings forming the junction. Note that we have kept 
the gauge choice (2.6) unperturbed by using worldsheet reparametrization invariance. We 
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then calculate the Nambu-Goto action for this ansatz and we expand it in powers of the 
fluctuations. The zeroth-order term gives just the classical action and the first-order 
vanishes thanks to the classical equations of motion.^ The resulting expansion for the 
quadratic fluctuations for the i = 1,2 string is written as [13] 



S. 



1 

"2^ 



dtdu 



fx 

2F- 



6x? - 



1/2 



gfy 



2E 



,3/2 



/2 



^gfy 
h 



2F 



1/2 



(3.2) 



+ 



fe 



y2 



hfeF, 



1/2 



2F. 



1/2' 



4F. 



1/2 



where we have used, the condition (2.9) and all functions and their ^-derivatives are 
again evaluated at = ^^q- Writing down the Euler-Lagrange equations for this action, 
using independence of the various functions from t and setting 



5a;f (t, u) = 6xf{u)e 



—iuit 



1,2,3, 



(3.3) 



we find that for the strings 1 and 2 the linearized equations for the three types of fluctu- 
ations read 



fx d 
du \ p^/^ du 

gfy d 



du \ p^/"^ du 



d 



fe d 



du \ p^/"^ du 



+ a;- 



MxF, 



1/2' 



gfy 



Sxi — , 



h 



,1/2 



Syi = o , 



Mof; 



1/2 



gfy 



2F, 



1/2 



f pV2 
o2 _ JyO^i o2 f 

SJy 



For the straight string the action for the quadratic fluctuations is 



:(3) 



47r/ 



dtdu 



^/g 



fe 



feh 



\fg \fg 



■0. 



h 



Sy 



fxh 



\fg \fgfy 



5x1 



wl + 



^/gfy 4^ 



(3.4) 



(3.5) 



careful treatment of the first-order terms gives rise to boundary terms. Demanding that they 
vanish turns out to be equivalent to the zero-force condition at the junction point. This can be shown 
quite straightforward for junctions in a flat spacetime before a gauge choice is made and reproduces the 
conditions found, for instance, in [10]. However, in curved spaces after the gauge choice u = ct is made, 
one obtains straightforwardly only the ^/-component of the zero-force condition. Switching to the gauge 
y = a one obtains the u-component as well. 



10 



from which we obtain the equations 



/ Vdfy 



du 



^ f fy £ '\ , 2 ^ 



A A^^'^ I ^.2 ^/e ^ 1^5^3 

Note that for the straight string the equations can be obtained by formally letting —>■ g 
and fyo ^ in (3.4). 

Hence determining the stability of the string configurations of interest has been reduced 
to a eigenvalue problem of the general Sturm-Liouville type 



A. 
du 



r{u; uq) > = (jj^q{u; uo)^{u) , u^in ^ uq ^ u < oo , (3.7) 



where the functions p{u]Uq), q{u]Uo) and r{u;UQ) are read off from (3.4) and (3.6) and 
depend parametrically on Ui through the function Fi in (2.12), which in turn is determined 
by (2.16) in terms of the junction point Uq and the NS and RR charges of the strings. 
Our aim is to find the range of values of uq for which cu"^ is negative. In fact the extremal 
such value of uq (the minimum as it turns out) will be decided by determining the zero 
mode a; = 0, which is an easier problem. Note also that, although it may seem so, the 
equations above for the three strings and for given type of fluctuations are not decoupled. 
The matching conditions at the junction point couple them in an essential manner. This 
equation can not be solved analytically in general, so that sometimes it is convenient to 
transform the Sturm-Liouville into a Schrodinger equation and use known approximating 
analytical methods that boil down to simple numerical problems. 



3.2 Boundary and matching conditions 

To fully specify our eigenvalue problem, we must impose boundary conditions on the 
fluctuations at the UV limit — > oo, at the IR limit u — Uj^m, as well as matching 
conditions at the junction point at u — Uq. The boundary condition at the UV are 
chosen to be 

$(it) = , as It ^ oo , (3.8) 

where $ is any of the fluctuations for the strings 1 and 2. These represent the fact that 
we keep the dyons at fixed points at the boundary. In the far IR at u — Umin we simply 



11 



demand, for the string 3 that extends also in there, finiteness of the solution and its 
M-derivative. Otherwise, all perturbations might be driven away from small values. 

We demand that the 5x fluctuations as well as their variations in the quadratic actions 
should be equal at the junction point u — Uq, so that the latter does not break. Moreover, 
we demand that the total boundary term resulting in deriving the classical equation of 
motions from the quadratic actions vanishes. These two requirements give rise to 

5xi — 5x2 — Sxs , at u — uo , 

Tp^q cos 9i5x'i + Tp/^q/ 008626x2 — Tm,nSx'^ = , at U = Uq , (3.9) 
where we have used (2.16) to simplify the second condition. 

Identical conditions with those in (3.9) hold for the angular fluctuations 59 as well 

69i = 602 = 603 , ai u = Uq , 

Tp^q COS ei5e[ + Tp>^q, COS 92592 - T^,„5^3 = , at u^uq . (3.10) 

The appropriate boundary and matching conditions for the 5y fluctuations should be 
found flrst in a coordinate system in which the classical solution does not change. This 
is simply given by [13] 

u — u -\- 5u{t,u) , 5u(t,u) — — ^T7~r ^ (3-11) 

VclW 

and is easily checked that the classical solution is not perturbed at all. This redeflnition 
does not affect the 5x- and (5 ^-fluctuations since they have trivial classical support and 
we keep only linear, in fluctuations, terms. At the junction point the 5u fluctuations 
are continuous and this leads to a discontinuity for the 5y fluctuations due to (3.11) 
for the strings 1 and 2. As before, we demand that the total boundary term resulting 
in deriving the classical equation of motions from the quadratic actions 82^ vanishes, 
keeping in mind that the variations of the 5?/j's in their quadratic actions should be equal 
at the junction point. With the use of (2.16) we have the following conditions for the 5y 
fluctuations 

5yi cot 9i + 5y2 cot 92 = , at u = uq , 

Tp^q cos^ 9i5y[ + Tp,^q, cos^ ^2(5^2 - ^m.nf^ya = > at u^uq . (3.12) 
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3.3 Zero modes 



As we have noted we can obtain the boundary of stabihty in parametric space by studying 
the zero-mode problem. Hence, we examine closely the three distinct types of fluctua- 
tions. 



3.3.1 Transverse zero modes 

For the transverse fluctuations the zero mode solution obeying (3.8) for the strings 1 and 
2 is 

Sx, = aMu), ^u)=r^J-^^, t = l,2, (3.13) 

Ju Jx ^ Jy Jyi 

where the a^'s are multiplicative constants. The zero mode solution for the straight string 
is 

5x3 — Ci3 du^ — h const. (3-14) 

Ju Jx 

When we specialize this expression in our examples we will see that either Sxs and/or 
Sx'^ diverges when u — > Wmin, so we choose 03 = in order for the solution and its 
derivative to remain finite. Thus Sx^ — const, and the zero mode solution becomes 
irrelevant for the stability analysis. Prom (3.9) and the zero force condition, one finds 
a hnear homogeneous system of algebraic equations for ai and 02, whose determinant 
should vanish for non-zero solutions to exist. In this way we find the following condition 

sin ^1/1 (mo) + sin ^2/2(^0) = . (3.15) 

Since Ii{u) > this condition cannot be satisfied and therefore the transverse zero modes 
do not exist, rendering the corresponding fiuctuating modes as stable. Our finding is 
similar to the conclusion that the transverse fiuctuations of the single string fluctuation 
corresponding to the potential of heavy quark- ant iquark pairs are also stable [13]. 



3.3.2 Longitudinal zero modes 

We next turn to the longitudinal fiuctuations. The zero mode solution obeying (3.8) for 
the strings 1 and 2 is 

Syi = biJi{u), J^(u)= I du _ ^ , i = l,2, (3.16) 



13 



where the ftj's are multiphcative constants. The zero mode solution for the straight string 
is 

Sys^h / du^ + const. (3.17) 

Ju Jy 

As before, when we speciahze this expression in our examples we will see that either 5yz 
and/or diverges when u — > Wmin, so we choose 63 = in order for the solution and its 
derivative to remain finite. Therefore the zero mode fluctuation becomes irrelevant for 
the stabihty analysis. Similarly to before, from (3.12) and the zero force condition one 
flnds the following relation 

cos9iJi{uq) — cos92J2{uq) . (3.18) 

This is an equation for Uq in terms of the string coupling and the strings' charges, which 
can be solved numerically. We will denote its solution, whenever it exists, by uqc- In 
general we will have expansions of the form 



n=l 

00 

n=l 

In the most important case of a quark-monopole interaction, S'-duality invariance of the 
result dictates that UQc{gs) — u^dgj^) and therefore the two perturbative expansions 
above contain the same information. Namely, Uq} — u^f^ and all coefficients are equal, 
i.e. XLfi — Vjji . Actually, then the critical value uqc does not change significantly from 
its leading order values u'q} , the reason being that the small and large Qs results actually 
coincide and therefore there is not much room for big changes in between. We show 
below that the leading term uf} can be computed by (3.24) which is an approximative 
version of (3.18). In the examples below the leading order result for uqc from (3.24) and 
that computed for Qs — ^ using (3.18) differ by two to ten percent. 

A natural question is whether or not the value of uqc coincides with that corresponding to 
the maximum value of L{uq), labelled by i^om if that exists, below which the energetically 
less favorable branch in the E — L diagram starts (a typical case is presented in fig. 
2). This might have been expected since in [13, 14] it was shown quite generally that 
an instability of longitudinal fluctuations of a single string used to compute the heavy 
quark- ant iquark potential, occurs precisely at the solutions of the equation L'{uq) = 
(which is identical to E'{uq) = 0). In our case following similar manipulations as in [13] 
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(a) (b) 

Figure 2: Plots of L{uo) and E{L) for cases where there is an instabihty of the longitudinal 
modes. This covers the case of non-extremal D3-brane, multicenter £)3-branes on a sphere 
for the ^0 = trajectory and the Rindler space. The various types of lines correspond 
to stable (sohd dark), metastable (sohd gray) and unstable (dashed gray) configurations, 
with the stability determined by the analysis of section 5. 

and using (2.16) at the equilibrium point where the 9iS are constant, in order to find the 
derivatives of Ui with respect to Uq, we obtain 

f 

L'{uo) = (sin^iXi + sin^si^s) , 

V Jyo 



E'(uo) = ^f'sme,Tp,,(smeiK^+sme2K2), (3.20) 
which are indeed proportional to each other and where we have defined^ 



Ki^ I dud^\y^\- i = l,2. (3.21) 

J Uo 



The value of uqc for which (3.18) is satisfied is different from uq — uom for which L'{uo) — 
E'{uq) — 0. In fact, in the examples we will explicitly work out, it occurs at smaller 
values, which implies that part of the upper branch of the solution is still perturbatively 
stable although energetically not favorable. However, the disconnected configuration 
becomes energetically favorable for all positive values of the energy. Thus part of the 
upper branch which is perturbatively stable is in fact metastable. The fact that Uqc 7^ 
can be demonstrated quite generally in the small coupling constant limit for the most 



^In cases, such as in Rindler space, where the upper Umit of integration is finite, i.e. «max = A, then 
there appears an additional term in the expressions below given by —^/Fi/f'y computed at it = A and 
also we replace infinity by A in the upper limit of integration. 
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significant case, namely for tlie quark-monopolc system. We may sliow tliat 

L'iuo) = Tdu (^] \ + 0{g,) , (3.22) 

y/jyo Juo \ ■'y J ~ ^y^' 

with an analogous expression for E'{uq). From this we compute the leading term u^q^^ in 
an expansion of the values Mom for which the maximum energy and length occur analogous 
to (3.19). Using (2.22), (2.23) and the general identity 

Ji{uo) = 2K,{uo) + 2 ^ , (3.23) 

proved by partial integration, we may show that the zero- mode condition (3.18) becomes 

2^- rdu^^Oig;), (3.24) 

JyO J uo Jy 

from which one computes the leading order value of Uqc in the expansion (3.19), namely 
Uq'J. Coming from different conditions we have that generically u'^ ^ Uq'J and this 
explicitly demonstrates that the value for Uom{gs) found by requiring maximum length 
and energy is different than the value of Uoc{gs) found by requiring the existence of a 
zero mode for the longitudinal fluctuations which signals perturbative instability of the 
system. 

We note that the condition (3.24) can be directly obtained from the boundary conditions 
(3.12) specialized to the zero mode solution for the longitudinal fluctuations. In the 
derivation one should be careful and absorb a factor of l/gg into the overall amplitude 
of the string 2. In this way it becomes apparent that although in the small coupling 
limit the string 2 becomes stiff and approximately straight, its ultra small in strength 
fluctuations couple to those of the string 1. Similar considerations can be given for the 
higher modes as well. The above comments explain also why although in the limit gg ^ 
(or gs oo) the turning point Ui uq (or U2 — >■ Uq) and therefore the corresponding 
equation for the fluctuations (3.4) becomes the same as the equation that appears in 
the case of the single string for the quark- antiquark system [13], the conditions for the 
existence of the zero mode are different in the two cases. They are so due to the different 
boundary conditions. 

Finally, we comment on whether or not the heavy dyon-dyon potentials as computed 
within supergravity using the AdS / CFT correspondence violate some general principles 
regarding the strength and sign of the corresponding force. For instance, for the quark- 
antiquark potential it turns out that the force should always be attractive with strength 
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an ever non-increasing function of the separation distance [17]. This is a concavity 
condition for the quark- ant iquark potential and imphes that the upper branch in the 
E—L plot should be disregarded, in full agreement with the results of [13, 14] employing 
perturbative stability. In our case, using (3.20), we can show that 

dE _ Tp^q sin Oi 



dL 2n 
d'^E _ Tp^qsin6i fyQ 



fyO , 

(3.25) 



which can be written in a symmetric way under an exchange of the strings 1 and 2 
using the y-component of the non-force condition. The expressions (3.25) for the dyon- 
dyon potential are analogous to those found in [16] for the quark-antiquark potential. 
Hence, the force is always attractive, but for the upper branch in the E—L plot in which 
L'{uq) > 0, has an increasing strength. On the other hand we have found that part 
of the upper branch is perturbatively stable, so that requiring concavity seems to be a 
stricter condition than perturbative stability. However, concavity is not a condition that 
gauge theory analysis of Wilson loops for heavy quark-monopole pairs requires. One 
can derive using reflection-positivity, that the energy of this pair is larger or equal than 
the average of the energies of the quark-antiquark and monopole-antimonopole systems.^ 
This condition is rather trivially satisfied in our case, provided we use the physical lower 
branches of the two pairs. 



3.3.3 Angular zero modes 

For the angular zero modes we cannot explicitly write down the solution for the zero 
mode due to the presence of the restoring force term in the corresponding Sturm-Liouville 

equation in the third lines of (3.4) and (3.6). Therefore, in this case we have to work 
out explicitly the result in the various examples we will present. As it was shown in 
previous work [13, 14] one can use approximate analytic methods in which the Schrodinger 
description of the differential equation (3.7) plays an important role. 



4 Examples of classical string junction solutions 

In this section, we study first the behavior of the dyon-dyon potentials emerging in 
Wilson-loop calculations for non-extremal and multicenter D3-brane backgrounds. Since 
■^We thank C. Bachas for providing this information. 
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we will work with the Nambu-Goto action, we need mention in the expressions below 
only the metric and not the self-dual five-form which is the only other non-trivial field 
present in our backgrounds. At the end we work out the example of the Rindler space 
which captures the behavior of black holes in various dimensions near the horizon. 



4.1 Non-extremal D3-branes 

We start by considering a background describing a stack of N non-extremal D3-branes 
[18]. In the field-theory limit the metric reads 

+ ^'(^7^ ^^' + ^^5) , (4.1) 



where the horizon is located at u = fi and the corresponding Hawking temperature is 
T — This metric is just the direct product of AdSg-Schwarzschild with and 



it is dual to jV' = 4 SYM at finite temperature. For the calculations that follow, it 
is convenient to switch to dimensionless variables by rescaling all quantities using the 
parameter /i. Setting u ^ /lu and uq — > /iuq and introducing dimensionless length and 
energy parameters by 

L^-L, E^^E, (4.2) 

/J, ZTT 

we see that all dependence on /i drops out so that we may set ^ 1 in what follows. The 
functions in (2.2) depend only on u (refiecting the fact that all values of 9 are equivalent) 
and are given by 

g{u) = 1 , fy{u) = (u^ - l)/i?^ , fM = {u' - 1)1 R' , 
4 1 4 

— 1 



= , h{u) = . (4.3) 



In this case u^^^^ — Wj-oot = 1) coinciding with the location of the horizon. For the dyon- 
dyon potential we find that the separation length is given in terms of the junction point 

Mo as 



L ^ Pi' 




du 



J^i — TT^ -f^l 



Zul H4'2'2'4'm4'«4 / ' 
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whereas the energy is given by 

E — Tr,^£■^ + 



= / du[. -^—,-l\-{u,-l) (4.5) 



1,2 



In the above Fi(a, 61, 62, c, Zi, ^2) is the Appell hypergeometric function and Uq'^ 1. 
Prom (2.16) one finds that the turning points of the two strings are given by 



Ui 



{ul-{ul-l)cos^eif/\ i = l,2, (4.6) 



where the angles 9i are given by (2.18), in terms of the strings' charges. From (3.20) we 
have in this case 

L'{uo) = i^' , ° (sin ^1X1 + sin 92K2) , 

K = -^F (I J_ <] _ J_f, /3 1 1 7 J_ 

' 28x^5 ' V4'2'2' 4 'u^'u^y' 12ul '\rm'ut'utJ ' ^ ' 

The function L{uq) has a single global maximum for any value of the string coupling. Its 
location depends on the string coupling and we obtain the following expansion for 
the quark-monopole system 

uom ^ 1.177 - 0.037^, + 0{gl) . (4.8) 

For L > -Linax , only the disconnected solution exists. For L < Lmaxj Eq. (4.4) has 
two solutions for uq, corresponding to a short and a long string and, accordingly, E is 
a double-valued function of L. The behavior described above is shown in the plots of 
Fig. 2 and we see that it is very similar to the behavior of the L{uo) and E — E{L) for 
the heavy quark-antiquark system [19]. However, as we have pointed out and we will see 
later in detail the upper branch, although energetically less favorable than the lower one, 
it is not unstable in its entirety. We note that the classical string junction in this case 
has also been examined in [20] where the behavior of fig. 2 was also noted. 



4.2 Multicenter D3-branes on a sphere 

Wc now proceed to the case of multicenter D3-branc distributions. These distributions 
[21, 22] constricted as extremal limit of rotating D3-brane solutions [23, 24] have been 
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used in several studies for the Coulomb branch of A/" = 4 SYM within the AdS/CFT 
correspondence, starting with the works of [16, 25]. Here, we will concentrate on the 
particularly interesting cases of uniform distributions of D3-branes on a three-sphere. 

The field-theory limit of the metric for N D3-branes uniformly distributed over a 3-sphere 
with radii Tq reads 

+ [{u^ - rl cos^ e)de'^ + cos^ Od^l + {u^ - rl) sin^ ^d^^ , (4.9) 

where ^ 

^ u^{u^ — rl cos^ 9) ' (^•-'■0) 
It is convenient to switch to dimensionless variables by rescaling all quantities using the 
parameter tq. Setting u — > TqU and Uq — > TqUq and introducing dimensionless length and 
energy parameters by 

L^-L, E^^E, (4.11) 
ro 27r 

we see that all dependence on ro drops out so that we may set ro — > 1 in what follows., 
we write the functions in (2.2) as 

g{u, 9) = ^^^^^ , fy{% e) = u'iu' - cos' 9)/R' , 

fa,{u,9) = u\u' -cos'^ 9)/ R"^ , (4.12) 

zt^ — cos^ 

f0{u, 9) ^ u' - cos^ 9 , h{u,9)^ ^— — — . 

U/ -L 

Therefore the conditions (2.9) are satisfied only for 9o = and 6*0 = 7r/2. We examine 
these two cases in turn. 

4.2.1 The trajectory corresponding to ^o = 

In this case Umin — liroot — 1 and the integrals for the dimensionless length and energy 
read 



R' I ujuj-l r '^'^ = + (1^2) 

= R'(^-^^[U{u,,k[',h)-F{u^,h)]+{1^2)^ (4.13) 
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and 



Si = r 

J Un 



du 



{u^ - ui){u^ + uj - 1) 



- {uo - 1) , 



(4.14) 



+ 1 



ul-l 



+ 1, i = l,2, 



where F(z/, A;), E(z/, A;) and !!(//,«,/;;) are the incomplete eUiptic integrals of the first, 
second and third kind respectively, while K(/c),E(/c) and Il{a,k) are the corresponding 
complete ones and 



Ui 



Hi — sm 



Prom (2.16) one finds that 



-1 / ^Q-^i 



i/j — sm 



2^.^ - 1 



(4.15) 



Ui 



l + y/Ai 



Ai ^ 1 + Auo{uq - 1) sin^ Oi , ^ = l,2, 



(4.16) 



where the angles 9i are given by (2.18). The function L{uo) has a single global maximum 
depending on the string coupling gg. The situation is similar to the quark- antiquark 

potential in this background [16] and to the case of string junctions on black DS-branes 
in the present paper. The general behavior is shown in the plots of Fig. 2. 



4.2.2 The trajectory corresponding to 9q — t:/2 

In this case Wmin = 1, but Wroot = and the integrals for the dimensionless length and 
energy read 



du 



'uo U-\/{u^ — l)(li^ - uf) 

i?2 fn{u,,i,ki)-F{u,,ki) 



+ (1-2) 



V2 



Ul 



+ (1 - 2) , 



(4.17) 
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(a) (b) 

Figure 3: Plots of L{uo) and E{L) for cases where there might be angular instabilities. 
Note the appearance of a confining potential for cases where there are angular instabil- 
ities. This covers, among our examples, the case of multicenter £)3-branes on a sphere 
for the 9q — t^/2 trajectory. The meaning of the colors is as in fig. 2. 



and 



"'^ duu I 



f 



1 - 



uo 



duu 



1 ' 



(4.18) 



ul-l 



where now 



'4+1 
2^2 



, k' = , 



//j = sm 

From (2.16) one finds that 



1 I i ^l-u^ 
ul-l 



2m? 



— sm 



ul + uj 



(4.19) 



Ui = Uq 



\/ sin 9i , i = 1, 2 , 



(4.20) 



where 9i angles are given by (2.18). In this case L{uo) is a monotonously decreasing 
function which approaches infinity as Mq — > 1 and zero as — > oo and hence no maximal 
length exists, Eq. (4.17) has a single solution for Uq given the length L and consequently 
i5 is a single-valued function of L. This behavior is shown in Fig. 3. This confining 
behavior is similar to what was found for the quark- antiquark system in [16]. 
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4.3 Rindler Space 



Our next example is based on the Rindler space, which it is a portion of fiat Minkowski 
space. In this space an observer experiences the Unruh effect, perceiving the inertial 
vacuum as a state populated by a thermal distribution of particles at a temperature 
T = K/27r, where k is the surface gravity. Black hole solutions near the horizon behave 
like the Rindler space, so that the study of string in this background serves more general 
purposes. In [26, 14] the string configuration corresponding to an open string with its two 
endpoints located at the same radius u — A (instead of w — > oo) we studied. Moreover, we 
have proven that this problem is exactly equivalent to the classical-mechanical problem of 
the shape of a soap film suspended between two circular rings and we have also performed 
a small fiuctuations analysis around this classical configuration. 

The metric for Rindler space has the form 

ds^ ^ -i^^u^dt^ + dy'^ + du"^ + . . . , (4.21) 

where u is the radial direction with the Rindler horizon corresponding to u — and y 
is a generic spatial direction. It is natural to extend the computation of [14] to the case 
of sting junctions with the end points of strings 1 and 2 at u — A and the third straight 
string stretching to u — 0. For convenience we pass to dimensionless units 

Au , uo^ Auo , AL , E ^ — -E (4.22) 

and then we use the formalism of section 2 which readily applies to this setup. The 
length and the energy of the string are determined by Eqs. (2.13) and (2.14) (without 
the subtraction term for the energy and the upper limit in the integration changed to 1) 
which for our case read 

C du ^ + (1 ^ 2) = fcosh-^ — - cosh-^ + (1 ^ 2), (4.23) 

J no \/v? -Ul \ Ui UiJ 

and 

E = Tp^q^i + Tp>^g>S2 + T„j^„7rMo , (4.24) 
= / f^^^ . = (jl-u] + uf cosh"^ — - uoJul-ul - ul cosh~^ — ) . 

Juo V ^ -^i \ Ui 1 UiJ 

Prom (2.16) one finds that 

Ui — u^smOi , i — 1,2 , (4-25) 
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where the angles 9i are given by (2.18). The function L{uq) has a single global maximum. 
Its location, depends on the string coupling Qg. For small values of the latter we get the 
expansion 

- 0.552 - 0.450y, + 0.47%^ - O.lOO^f + 0{gf) . (4.26) 
The situation is similar to that depicted in the plots of fig. 2. 

5 Examples of stability analysis 

In this section we apply the stability analysis developed in section 3 to the string solutions 
reviewed in section 4. We find that for the non-extremal D3-brane, the sphere with = 
and the Rindler space we have a longitudinal instability corresponding to the upper 
branch of the energy curve, for the sphere with = ti'/S we have angular instabilities 
towards the IR, even though the potential has a single branch. 

5.1 The conformal case 

For completeness we consider first the conformal case, corresponding to the /i — or 
ro ^ limit of any of the above solutions. The classical string junction solution has been 
constructed in [12]. Since there is no ^-dependence, the 59 fluctuations are equivalent to 
the transverse ^^-fluctuations and therefore stable. On the other hand for the longitudinal 
fluctuations we have to solve (3.18). From (2.16) one finds that 

Ui = KoVsin^i , i = l,2, (5.1) 

where the angles Oi are given by (2.18). Moreover we compute 

3^"^ (I'll- 1)- (^-2) 

Using the above, we may check that (3.18) has no solution for any value of string coupling 
Qsi thus we verify that indeed string junctions are stable in the conformal case. 

5.2 Non-extremal D3-branes 

We proceed next to the case of non-extremal D3-branes, where we recall that the potential 
energy is a double- valued function of the separation length. As before, since there is no 
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^-dependence, the 66 fluctuations are equivalent to the dx-fluctuations and stable. For 
the longitudinal fluctuations we have to solve (3.18), where 



Equation (3.18) has a solution for any value of string coupling Qs- This value Uqc, is 
smaller than the value uom in which the maximum of the length occurs. In the limit 
of small coupling constant for the quark-monopole system we have shown that we can 
approximate (3.18) by the condition (3.24) which in our case becomes 

Its solution gives the leading order result for the critical value in which instability occurs 

uoc ^ 1.117 - 0.166^, + 0{g'',) , (5.5) 

where we have included the flrst correction as well. This expansion is different than 
the expansion of Uom in (4-8). Thus the whole upper branch is not entirely unstable, 
a behavior not expected from energetic considerations. Let's compare the maximum 
length of the screened interaction for the cases of the quark- antiquark and the quark- 
monopole interactions. As gs varies from small to large values the maximum length is 
- (0.425±0.005)ii:2, whereas for the quark-antiquark pair L^^ ~ 0.8697?^ ~ 2L^. 
We see that the quark-monopole pair is screened twice as strong as the quark-antiquark 
pair. This is a general feature in all of the examples we encounter in this paper. 

5.3 Multicenter D3-branes on a sphere 

The results of the stability analysis for the two orientations for the sphere distribution 
are presented below. 

5.3.1 The trajectory corresponding to = 

Angular fluctuations are stable due to the fact that, as in [13], we may transform, for 
each string separately, the problem into a Schrodinger equation with a positive potential 
for all values of Ui. For the longitudinal fluctuations we have to solve (3.18), where 

M'^o) = I du — TTT^ ' i = 1,2 . (5.6) 



uo 



(^2(^2 - 1) - ^2(^2 _ i))3/2 ' 
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This has a solution for all vahics of the string coupling. In the limit of small string 
coupling for the quark-monopole system the condition (3.24) reads explicitly 



2ul Uq . f Uq + I 



^Inr^ =0{9,), (5.7) 



2ul-l 2 V^*o - 1, 
giving the leading order result for the critical value in which instability occurs 

uoc 1.08A + 0{gs) . (5.8) 

Thus, part of the upper branch of the solution is perturbatively stable as in the case of 
the string junction based on the non-extremal D3-brane solution. As before, we compare 
the maximum length of the screened interaction for the cases of the quark- antiquark and 
the quark-monopole interactions. As Qs varies from small to large values the maximum 
length in the quark-monopole interaction is ~ (0.425 ± 0.050)i?^, whereas for the 
quark-antiquark pair L^^ ~ 1.0027?^ ~ 2.3L^^. 

5.3.2 The trajectory corresponding to Oq = tt/2 

Longitudinal fluctuations are stable since an explicit check of (3.18) shows that it has no 
solution for any value of string coupling gg. Alternatively, we may transform the problem 

into a Schrodinger problem as in [13] with an everywhere positive potential. For the 
angular fluctuations the potential is just constant 

Ve^-l, i = l,2. (5.9) 

The change of variables for the two string is 

Zi{u) = F{ui, k) , i = 1, 2 , (5.10) 

UiV2 

where Ui and ki are given by 



= l^/ .^ + .gsingJ ' ^-^V sin ' ^'-''^ 

The zero mode of the Schrodinger equation for the i = 1, 2 strings read 

= Cisinzi , Zi e [0,0] (as u e (oo,mo]) , (i = ^i(«o) , (5-12) 

where q are multiplicative constants. The straight string 3 plays a role in the bound- 
ary conditions, the reason being that the mass term is present in the equation for the 
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angular fluctuations even for the zero mode. For the straight string we have two linear 
independent solutions sin^ and cos 2;, where the appropriate change of variable is 

z = - - sin-i - . (5.13) 

We examine in detail the case where only sin 2; is used and comment on the general 
linear combination later. Prom the matching conditions in (3.10) one finds the following 
condition 



sin 9i sin G cos C2 + sin 62 sin G cos Ci — ^^'^l^i"*" sin Ci sin C2 = . (5.14) 



For the quark- monopole case we may use that + ^2 = 7r/2, to simplify it a bit. Fur- 
thermore, in the small string coupling limit (5.14), becomes 

which is a more explicit and easy to handle transcendental equation. Hence, the solution 
for the leading order critical values for which instability of the angular perturbations 
occur is 

uoc 1.378 + 0{gs) . (5.16) 

The confining behavior turns out to be unstable since (5.14) has a solution for any value of 
the string coupling Qs- As Qs varies from small to large values, the maximum length of the 
quark-monopole interaction is ~ (0.560±0.040)i?^, whereas for the quark- antiquark 
pair Lf^ ~ 1.700i?2 ~ 3L^^. 

Let' s mention that, had we used the general solution of the Scrodinger equation, that is 
cos(2; + (p), we would have obtained similar results for any phase with the exception of 
(y9 = for which case no instability occurs. Hence the above considerations and results 
are quite generic. 

Finally let's note that, in order keep the discussion at a reasonable length, we refrained 
from explicitly presenting the details for the other important uniform distribution of 
D3-branes, namely, that on a disc. We mention the end result: For Oq — 7r/2 the solu- 
tion is stable against all types of small fluctuations, while for ^0 = it is stable against 
longitudinal fluctuations and unstable against angular ones after a certain separation 
length. This behavior reveals the expected screening behavior which moreover is practi- 
cally independent of the orientation of the string and similar to that for the case of the 
quark- antiquark potential studied in full detail in [13]. 
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5.4 Rindler Space 



Finally, we consider the case of the Rindler space. Fluctuations of the single string 
solution corresponding to the heavy quark- antiquark potential have been considered in 
[14, 26]. Turning to our case, we see from the form of the metric only the longitudinal 
fluctuations are relevant for our discussion, so that we have to solve (3.18), where 

Jiiun) = — I cosh^^ cosh~^ — , H " I . (5.17) 

Equation (3.18) has a solution for all values of the string coupling Qs and we have the 
same behavior as in 6y fluctuations in the non-extremal D3-branes. For small values of 
the string coupling latter we get the expansion 

uqc ~ 0.368 - 0.214^, + 0.218^,^ - 0.058^^ + 0{g^^) . (5.18) 

We see that it is different from the expansion of uom in (4.26) corresponding to the 
maximum of the length and the energy. Hence, the behavior is similar to that represented 
in fig. 3. 

6 Discussion 

In this paper we have constructed string junctions on general backgrounds and used them 
to compute the interaction energy of heavy dyons at strong coupling within the AdS/CFT 
correspondence. The most important case is that of the quark-monopole pair for which 
we paid particular attention. Then wc turned into the pcrturbative stability analysis 
of these solutions which is important to distinguish the physically relevant regions in 
which the potentials can be trusted. We formulated and further investigated general 
conditions for the existence of instabilities. Starting with the conformal case we indeed 
verified that there are no instabilities at the conformal point of jV = 4 SYM, as expected. 
Nevertheless, we emphasize here that in the /3-deformed conformal solution of [27], as 
well as for all backgrounds that asymptote this, there is an instability for certain angular 
fluctuations. The reason is identical to the case of a single string for the quark- antiquark 
potential that was found in [14] (see section 4.3.1) and occurs for values of the real 
deformation parameter a (in the standard notation) larger than a certain critical value 
(we remark that the classical string junction of the /3-deformed background of the disc 
D3-brane distribution has been considered in [28] and should suffer from this instability) . 
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We found that part of the branches of the potential energy versus length that could be 
disregarded on energetic considerations can be perturbatively stable. The reason is traced 
to the fact that the matching conditions at the junction point allow for interaction of all 
strings in the junction that enhances its stability. This is unlike the findings in [13, 14] 
for the quark- antiquark potential according to which the perturbative instability renders 
energetically unfavorable branches as unstable. We pointed out that this result is not in 
conflict with general considerations on the quark-monopole potential. For instance, there 
is no concavity condition as is in the case of the quark- antiquark potential. We found 
that the confining behavior of the dyon-potential in the Coulomb branch of the J\f — 4: 
SYM with vacuum expectation values distributed uniformly on a sphere, is unstable. 
This is similar to the result of [13] for the seemingly confining behavior of the quark- 
antiquark potential for the same distribution. It is quite remarkable that demanding 
perturbative stability of the string configurations renders as physically unacceptable the 
confining behavior, which indeed is not expected from gauge theory considerations. 

Our general formalism was based on a diagonal class of metrics (2.1). It can be extended 
to supergravity backgrounds in which off-diagonal metric elements appear and in which 
additional background form fields couple in the Wess-Zumino part of the string actions 
in the junction, in analogy with [14] for the quark- antiquark case. In particular, this 
will cover, within the AdS/CFT correspondence, interactions of moving dyons in hot 
quark-gluon plasmas. We expect that also in this case the associated energy will have 
two branches with the upper one being only partly unstable. The behavior found in 
[29], namely that the maximum length behaves for high enough velocity as L^ax — 
0.743(1 — is expected to persists with the overall numerical coefficient replaced 

by, roughly, its half. 

It will be interested to extend our analysis to junctions involving more than three strings 
where we believe similar phenomena with the triple-junctions considered here will be 
found. A more interesting generalization is to consider systems of multiple external 
particles at finite temperature and examine to what extend they prefer, under small per- 
turbations to split, with raising temperature, into smaller systems. This possibility was 
argued for, based on energy considerations, in [30] who examined in detail the interaction 
energy for the system of a quark, a monopole and a dyon. 

Finally, it has been proposed that F- and D-superstrings could have been produced at 
the early Universe, then expanded at a cosmic size and that their network properties 
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by forming junctions make them distinguishable [31] from gauge theory cosmic strings 
(see, for instance, [32]). The no- force condition for a stationary junction to form, gives 
rise to a locally flat space time [33], much-like the case of single cosmic strings [34]. 
Small fluctuations of cosmic string junctions on flat spacetime have been analyzed [35] 
in closed analogy with corresponding work on superstring junctions (see, for instance, 
[32]). In view of the above remarks it will be very interested to study the formation and 
stability of cosmic size string junctions in cosmological backgrounds, by adapting the 
results and techniques of the present paper, to such cases. 
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A An analog from classical mechanics 

Longitudinal string fluctuations exhibit a behavior which is not expected from the en- 
ergetics of the configurations. However, this behavior can be found in the classical me- 
chanics problem of determining the shape of a thin soap film stretched between two rings 
appropriately, modified to serve our purposes. We briefiy review first the solution of the 
unmodified classical problem which is a catenary and can be found in standard textbooks 
[36]. In the notation of the present paper and of the appendix of [13] the cylindrically 
symmetric solution reads 



in terms of the integration constant parameter uq. The potential energy and length of 
the soap film are given by the following expressions 
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r{z) — Uq cosh 



z 



L L 
2 2 ' 



(A.l) 



L — 2uo cosh ^ — , 

Uo 





(A.2) 
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Hence, E{L) is a double- valued function of L with the energetically favorable branch 
having lower energy for a given distance, corresponding to a shallow catenary, whereas 
the upper branch corresponds to a deep catenary. There is also the Goldschmidt solution, 
with two disconnected soaps in the two rings with energy E — 2tt. Thus we have the 
same qualitative behavior as in the classical solution of the non-extremal D3-branes, 
depicted in fig. 2. The value of uq at the turning point is uom — 0.552 leading to 
the maximal length of — 1.325. The Goldschmidt solution becomes energetically 
favorable at a higher value of Uq ~ 0.826, corresponding to a smaller value of L ~ 1.055. 
The perturbative small fluctuation analysis has been performed in [37] and in [13]. The 
small fluctuations stabihty analysis of this solution gives the expected result that the 
upper, energetically unfavorable, branch is perturbatively unstable. 

We would like to add an extra term to the quadratic fluctuations (see Eq.(A-l) in ap- 
pendix A of [13]) that, to leading order, would not affect the energy of the solutions and 
then check whether part of the upper branch becomes perturbatively stable due to this 
modification of the action. In order to fulfill the above requirement imagine that we first 
let the soap film assume its shape given by (A.l) and then we turn on the perturbation 
which we take it to be quadratic in the normal fluctuations of the soap film surface. The 
normal fluctuations decouple from the two tangential perturbations which moreover are 
trivial [13]. This will not affect the classical solution and its energy, but will certainly 
modify the stability analysis. There are several such terms that we may add. Here we 
choose to present the two most natural ones. 

In the first we add a term to the potential that acts only at the minimum of the catenary 
Bit z — 



where the factor uq is introduced for convenience. For the normal perturbations, we 
define u — z/uq, we separate variables according to 



k>0 



(A.3) 



d7]{t,u,(j)) = $(M)e-*™e*'^'^ 



m = 0,±l,±2,... , 



(A.4) 



and eventually we end up with the Sturm- Liouville equation 




(A.5) 



subject to the following boundary conditions 




(A.6) 
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representing fixed endpoints. Equation (A. 5) can not be easily solved for non-vanishing 
values of tu, but we can obtain useful information for our stability problem by considering 
the zero-mode problem having uj = 0. The transformation x = tanhw turns (A. 5) to 
an associated Legendre equation with the general solution given by a linear combination 
of P{"(tanhH) (they vanish identically unless m — 0, ±1) and Q'^{ta,nh.u), with different 
coefficients to the regions left and right of 2; = 0. However, we know from the work 
in [13] that in the absence of the extra term (for k — 0), only for m = a zero mode 
exists for the value Uq ~ 0.552 we mentioned above. Since the modification is by an 
repulsive (5-function term the possibility of having zero mode solutions in the presence 
of it for \m\ ^ 1, is definitely excluded. Hence we concentrate to the case with m = 0. 
The continuity of the solution at z = 0, the discontinuity of its derivative, as read off 
from (A. 5), and the boundary conditions (A. 6), give a linear homogeneous system for 
the four independent coefficients which in order to have non-trivial solution should have 
a vanishing determinant resulting to the condition 



KPi{a) = 2Qi(a) , =^ -|- 2 = aln ^ | ^ J , a = \J 1 — Uq . (A.7) 

This equation has a solution for all k ^ 0, which varies monotonously from uqc — 0.552 

to 0, as K G [0, 00). Importantly, the product -Uqc/^ and therefore the term (A. 3) we have 
added, remain finite. For small and large values of k we have the behaviors 

uqc ^ 0.552 - 0.160k + 0{k^) (A.8) 

and 

i^oc = 2e-'^/2-i ^ _ (A.9) 

An alternative to (A. 3) potential term is 

V^^-^{5^)\ k^O, (A.10) 

where the factor Uq is, as before, introduced for convenience. In the Sturm-Liouville 
equation (A. 5) we replace k6{u) u^k. The same transformation x = tanh-u gives an 
associated Legendre equation with the general solution given by a linear combination of 
P["(tanhM) and Q'^{tanhu), where m — \/w? -|- u^k. For identical reasons as before we 
concentrate to the case with m = and in fact refining the argument we should have 
that m < 1 in order to find a normalizable zero mode. Imposing the boundary condition 
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(A. 6) and demanding that the determinant of the resulting Unear homogeneous system 
for two independent coefficients has non-triviai solutions, we obtain the condition 

pf^{a)Qf^{-a) = P/^{-a)Qf^{a) , a = sjl-ul . (A.ll) 

This can only be solved numerically and as before the product uock remains finite. For 
small and large values of k we have the behaviors 

Moc^ 0.552 - 0.102A;^ + C(A;^) (A.12) 

and 

u^c^\-^ + 0{l/k'') . (A.13) 

In both of the above examples we see that the critical value Uqc can get arbitrarily small, 
which implies that the entire upper, energetically less favorable branch corresponding to 
the deep catenary, can become perturbatively stable. 
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